NON-UNIQUENESS FOR NON-NEGATIVE SOLUTIONS OF 
PARABOLIC STOCHASTIC PARTIAL DIFFERENTIAL 
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Dedicated to Don Burkholder. 

Abstract. Pathwise non-uniqueness is established for non-negative so- 
lutions of the parabolic stochastic pdc 

— = —X + X?W + ip, X = 

at 2 

where W is a white noise, ip > is smooth, compactly supported and 
non-trivial, and < p < 1/2. We further show that any solution spends 
positive time at the function. 



1. Introduction 

Let a : R — > R be p-Holdcr continuous (so \a{x) — <r(y)\ < K\x — y\ p ), 
let ip £ Cc(K) (the space of C 1 functions on M. with compact support), and 
consider the parabolic stochastic partial differential equation 

(1.1) ^(t, x) = jX(t, x) + a(X(t, x))W(t, x) + 4k 

Here W is a space-time white noise on R + xR. If a is Lipschitz continu- 
ous, pathwise uniqueness of solutions to (jl.ll) is classical (see, e.g., jWa!86j ) . 
Particular cases of ()1.1|) for non-Lipschitz a arise in equations modeling pop- 
ulations undergoing migration (leading to the Laplacian) and critical repro- 
duction or resampling (leading to the white noise term). For example if 
<j{X) = \fx and X > 0, we have the equation for the density of one- 
dimensional super-Brownian motion with immigration ip (see Section III. 4 of 
[POT] ). If a(X) = y/X(l - X), ip = and X e [0, 1] we get the equation for 
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the density of the stepping stone model on the line |Shi88j . In both cases path- 
wise uniqueness of solutions remains open while uniqueness in law is obtained 
by (different) duality arguments (see the above references). The duality argu- 
ments are highly non-robust and fail, for example if <r(x,X) = y/ f(x, X)X, 
which models a critically branching population with branching rate at site 
x in state X is f(x,X). This is one reason that there is interest in proving 
pathwise uniqueness in under Holder continuous conditions on a, cor- 

responding to the classical results of |YW71j for one-dimensional SDE's with 
Holder 1/2-continuous diffusion coefficients. 

In [MP10] pathwise uniqueness for (jl.ljl is proved if p > 3/4 and in 
jMMPll] pathwise uniqueness and uniqueness in law are shown to fail in (jl.ll) 
when a(X) = \X\ P for 1/2 < p < 3/4. Here a non-zero solution to (II. ip is 
constructed for zero initial conditions and the signed nature of the solution is 
critical. In the examples cited above the solutions of interest are non-negative 
and so it is natural to ask whether the results in jMPIO] can be improved 
if there is only one point (say u = 0) where u(u) fails to be Lipschitz, and 
we are only interested in non- negative solutions. Finding weaker conditions 
which imply pathwise uniqueness of non-negative solutions in this setting is a 
topic of ongoing research. In this paper we give counterexamples to pathwise 
uniqueness of non-negative solutions in the admittedly easier setting where 
p < 1/2. Even here, however, we will find there are new issues which arise in 
our infinite dimensional setting. Our methods will also allow us to extend the 
nonuniqueness result in |MMPll] mentioned above to < p < 1/2. 

We assume W is a white noise on the filtered probability space (Q, J 7 , Tt,P), 
where Tt satisfies the usual hypotheses. This means Wt{4>) is an J^-Brownian 
motion with variance ||</>||| for each <fi € L 2 (R, dx) and Wt{4>\) and Wt{4>2) 
are independent if (</>i,</>2) = / 4>i(x)<p2(x)dx — 0. A stochastic process X : 
f2 x R + xl-)l which is Ft — prcvisiblc x Borel measurable will be called a 
solution to the stochastic heat equation (|1.1[) with initial condition Xq : R — >• M 
if for each G C~(R), 



(The existence of all the integrals is of course part of the definition.) It is 
convenient to use the space C rap (R) of rapidly decreasing continuous functions 
on K as a state space for our solutions. To describe this space, for / € C(M) 
(the continuous functions on R) let 




+ 




a(X(s,x))(j)(x)W{ds,dx) +£(</>,?/>) for all t > a.s. 



|/| A = sup e A N|/(:r)|, 
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and set 

Crap = {/ G C(R) : |/U < oo VA > 0}, 
Ctem = {/ G C(M) : |/U < oo VA < 0}. 
Equip C rap with the complete metric 

oo 

d(/,. 9 )=^2- fe (||/-.9|| fe Al), 
fe=l 

and Ct em is given the complete metric 

oo 

d tem (/,g) = ^2- fc (||/-.9||-i/ fe Al). 
fe=l 

Let C^j p be the subspace of non- negative functions in C rap , which is a Polish 
space. Our primary interest is in the smaller space C^ ap resulting in stronger 
non-uniqueness results. 

A C^-valued solution to (jl.lj) is a solution X such that t — > X(t, •) is in 
C(E + , C+j p ), the space of continuous C+ - valued paths for all u. In general 
if £ is a Polish space we give C (R+ , £7) the topology of uniform convergence 
on compact sets. 

The following result is proved just as in Theorem 2.5 of |Shi94j . 

Theorem 1. (Weak Existence of Solutions). Assume ip > and the p-Holder 
continuous function a satisfies u(0) = 0. If Xq G C^ ap , there exists a filtered 
space (H, J 7 , !F tl P) with a white noise W and a Cf ap -valued solution of p. II) . 

Proof. Our conditions on a imply the hypothesis on a in Theorem 2.5 of 
|Shi94j . however that reference assumes ip(x,X) satisfies ip(x,X) < c\X\. 
The proof, however, extends easily to our simpler setting of tp(x) > 0. □ 

Here is our main result on non- uniqueness. The proof is given in Section [3J 
Recall that ip G C£(R). 

Theorem 2. Consider dTTTJ with a(X) = \X\p for p G (0,1/2) and ip > 
J" ip(x) dx > 0. There is a filtered space (57, J 7 , J 7 *, P) carrying a white 
noise W and two C^ ap -valued solutions to with initial conditions Xq = 

Xq = such that P^X 1 ^ X 2 ) > 0. That is, pathwise uniqueness fails for 
non-negative solutions to ([l.ip for a, ip as above. 

Remarks. 1. The state of affairs in Theorem [5] for ip = but Xq non-zero 
remains unresolved. We expect the solutions to still be pathwise non-unique. 
The methods used to prove the above theorem do show pathwise uniqueness 
and uniqueness in law fail if ip = Xq = and we drop the non-negativity 
condition on solutions. Namely, one can construct a non-zero solution to 
the resulting equation. We will not prove this as stronger results (described 
above) will be shown in jMMPll] using different methods. 
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2. Uniqueness in law holds for non-negative solutions to for ip, a as 
above and general initial condition Xo £ C+, p but now with 1 > p > 1/2. 
This may be proved as in |My98| where the case ifj = is treated; for p = 1/2 
this is of course the well-known uniqueness of super-Brownian motion with 
immigration We do not know if uniqueness in law fails for p < 1/2. The 
presence of a drift will play an important role in the proof of Theorem [2] 

3. A key technique in this paper is to consider the total mass M t = (X t , 1), 
and then apply Theorem 4 which, given the Holder continuity of X(t,x), 
shows the brackets process [M] to be bounded below by the integral of a 
power of M . This in turn allows one to apply comparison arguments with one 
dimensional diffusions. 

In Section [4] below we prove that in the corresponding stochastic ordinary 
differential equation, although pathwise uniqueness again fails, uniqueness in 
law does hold. Of course the SDE is now one-dimensional so on one hand 
this is not surprising. On the other hand, the manner in which uniqueness 
in law holds is a bit surprising as the SDE picks out a particular boundary 
behaviour which has the solution spending positive time at (see Section @|. 
This leads naturally to the following property for all solutions to the SPDE 
in Theorem [2] 

Theorem 3. Assume a and ip are as in Theorem^ Let X be any C^ ap -valued 
solution to (jl.ljl with Xq = 0. Then 

t 

l(X(s, x) = Vx) ds > for all t > a.s. 
The proof will be given in Section [5] below. Let 

b = (Va i) > o. 

We note that the above result fails for p = 1/2 since in that case Y t = 4(X t , 1) 
is a Bessel squared process of parameter 46 satisfying an ordinary sde of the 
form 

dY t = 2^/YtdBt +4bdt. 

Such solutions spend zero time at (see for example, the analysis in Section 
V.48 of JEW].) 

Finally we state the non-uniqueness result which complements that in 
jMMPll] in the much easier regime of p < 1/2. The solutions here will 
be signed. 

Theorem 4. If < p < 1/2 there is a G ra p-valued solution X to 

BX AX 
(1.2) ?- { t,x) = —(t,x) + \X(t 7 xWW 7 X(0)=0, 

so that P(X ^ 0) > 0. In particular uniqueness in law and pathwise unique- 
ness fail in (|1.2p . 
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Although the construction in jMMPllj for 1/2 < p < 3/4 is more delicate, 
it is a bit awkward to extend the reasoning to p < 1/2 and so we prefer 
to present the result here. The proof of Theorem 0] is simpler than that of 
Theorem [2] in that we can focus on a single process rather than a pair of 
solutions. The two proofs are similar in that approximate solutions arc found 
by an excursion construction and the key ingredient required for the SPDE 
setting is Theorem [5] below. Hence we only give a brief sketch of the proof of 
Theorem 2] at the end of Section [3] 

2. A Real Analysis Lemma 

Theorem 5. If < a, /3 < 1 and C > 0, there is a constant K^f3, C) > 
such that if f : M — > R+ satisfies 

(2-1) \f(x)-f(y)\<C\x-yf, 

then 

J f a dx>I<^J fdx) [ 
Proof. First we use a scaling argument to reduce to the case 

f{x)dx = 1 

for which we would have to prove 

Indeed, if we take b > and let 

g(x) = b-tfQbx) 

then 

\g(x) ~ g(y)\ = b^\f(bx) - f(by)\ < C\x - yf 
by the conditions of Theorem [SJ Then setting y = bx, we get 

g(x)dx = [ b~ p f{bx)dx = b' [l3+1) { f(y)dy = 1 



(a/3+l)/(/3+l) 



provided 



b=(J f(y)dy 



So g satisfies the conditions of Theorem [5] with J g = 1, and if we could show 
that 

J g a {x)dx > K A 
it would follow, substituting for 6, that 

f a {x)dx = b afi+1 I g a (x)dx >K A [ I f(x)dx ' 
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as required. 

Now we concentrate on proving J f a > K4 assuming that J f = 1 and 
assuming the Holder condition (|2.1j) on /. Let M = s\rp x f(x), and note the 
conclusion is obvious if M = 00 so assume it is finite. If M < 1, then since 
< a < 1, we have 



f a (x)dx > J f{x)dx = 1. 

On the other hand, if A/ > I, then the Holder condition on / implies that 
/ > i on an interval / whose length is bounded below by a constant L > 
depending only on C, (3. So in this case, too, we conclude 

J f a (x)dx > ^ > |, 

and Theorem [5] is proved. □ 



3. Proof of Theorem O 

If V> S C£(R), V» > 0, 6 = / V'dx > and < p < 1/2, we want to construct 
distinct solutions X, Y to 

dX \X 

(3.1) — (t, x) = — (f, x) + x)YW(t, x) + ^(x), X>0,X = 0. 

Let denote the space of bounded C k functions on R with bounded j th order 
partials for all j < k, and set Cf, = C®. The standard Brownian semigroup is 
denoted by (Pt,t > 0) and p t {-) is the Brownian density. 

Here is an overview of the proof. We will proceed by constructing ap- 
proximate solutions (X s , Y £ ) to (|3.1[) and then let (X,Y) be an appropri- 
ate weak limit point of (X Sn ,Y £n ). These approximate solutions will satisfy 
X s > Y £ > and Y s > X s > 0, respectively, on alternating excursions away 
from by M = (X E , 1) V (Y 6 , 1). M will equal 2be at the effective start of 
each excursion. We then calculate an upper bound on the probability that 
M will hit 1 on a given excursion (see (|3.43p below) and a lower bound on 
D e = \(X e ,l) — (Y e ,l)\ hitting an appropriate xo £ (0,1) during each ex- 
cursion (see p.58|) below). Theorem [5] is used in the proof of the first bound 
(see p.36|) below). These bounds will then show there is positive probability 
(independent of e) of D e hitting xo before M hits 1. The result follows by 
taking weak limits as e n I 0. The use of Theorem [5] will mean the above 
upper bound is valid only up to a stopping time V£ which will be large with 
high probability. This necessitates a "padding out" of the above excursions 
after this stopping time, and this technical step unfortunately complicates the 
construction. 

Fix e > and define {X s , Y E ), D e = \X £ - Y £ \, the white noise W and 
a sequence of stopping times inductively on j as follows. Let Tq = 0, U £ = 
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T? + e, and assume Xf, e = Yf E = on {TJj- < oo}. Assuming {TJj < oo}, 
on [T^.C^.] define 

(3.2) r e (i,a;) ee 0, and D £ {t, •) = X e (t, ■) = 2 / P s V(0^ € C+ . 

Jo 



That is, 



9X £ A 



(3.3) __ = _x £ + 2^ forT|.<t<C/|.. 

Next let t -> (^7|. +t ,-D[r|. +t ) in C(R+,C*+ p ) 2 solve the following SPDE 

for 

(3-4) ^ = ^ + j> + (Y°YW, Y{j h ee 0, 

^ = ^ + [(F e + £» e f - (F s f]^,D £ as .(x) = 2^P s ^{x)ds. 

The existence of such a solution on some filtered space carrying a white noise 
follows as in Theorem 2.5 of [Shi94j . To be careful here one has to construct an 
appropriate conditional probability given Txj* . and so inductively construct 
our white noise along with (Y e , D e ). Set Xf = Y t e + D e t for < t < T$ j+1 , 
where 

T| j+1 = mf{t > W 2J : {XI I) = 0} (inf = oo), 
and also restrict the above definition of (Y £ , D € ) to [U^j, T , ; f 7 - +1 ] . Therefore on 

dY e A 

_ = -y* + v + (Y*yw, ys s , = o, 
9jf e a r e 

(3.5) — = -X F - + i, + (X e fW, Xfj^ (x)=2j P s ^(x)ds, 

D e = \X e - Y e \ = X e -Y e > 0, 
X s ,Y E arc continuous and C+^-valued. 
Note that X£, e = Y£ e = 0. The precise meaning of the above formulas 
for dY £ /dt and dX e /dt is that equality holds after multiplying by e 
and integrating over R and over any time interval in [t/J^ , 2^fj+i] ■ 

Now assume T% j+1 < oo and construct (X 6 ,Y E ) and D E = \X e - Y e \ = 
y e — X 6 as above but with the roles of X and Y reversed. This means that 



on [Ti J+1 ,Ui J+1 ] = [Ti j+1 ,Ti j+1 +e], 



(3.6) D £ (t, •) = Y e (t, •) = 2 / PM-) ds £ C+ p (R) and X £ {t, •) ee 0, 

Jo 

and so 

dY E A 

0.7) ^r-2 y£ + 2 ^ 
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and on [U$ j+1 ,T$ j+2 ], 

^ = ^x £ + ^ + (x°yw, x^ +i ee 0, X- > 0, 

dY E A ■ [ E 

(3.8) — = -F £ + V; + (Y'fW, F f f J . +i = 2 / P s ^(x) ds, F e > 0, 

l> £ = |x £ - r £ | = r £ - x £ > o 

X £ ,F e arc continuous and C+, p -valued. 
Here, as before, we have 

Ti j+2 =M{t>U! j+1 :(Y t s ,l)=0}. 

Clearly Xj,e a = ^xf. 2 = on {^2j'+2 < °°} anc ^ ^7 t 00 an d so our 
inductive construction of (X s , Y £ ) is complete. It is also clear from the con- 
struction that if Xj(t) = ^(V=+t)AT e an d similarly for Yj, then we may 
assume 

(3.9) P((X 2j ,Y 2j ,Ti j+1 - T 2 j) G 

= P((X ,y ,3i) e •) a.s. on {Tfj < oo}, 

and 

(3.10) P((Y 2j+1 ,X 2j+1 ,Ti j+2 - G -| Ft S]+1 ) 

= P((X ,Y ,T £ ) G •) a.s. on {T| J+1 < oo}. 

Define Je = \J*LAUt-i,T?}, 

oo „t 

A\{t,x) = ]T(-1) J / l[T;,u;](s)iKx)ds, 

j=0 Jo 

and A e 2 {t,x) = -A\{t,x). Combine (pO]) . (|3~5)l . (pT5|) and the fact that on 
p2j+i> we nave -^t = = Y^t e + ^ — V 1 to see that for a test function 

(3.11) (X t £ ,0) = (A\(t),4>) + f Uxi,^4>) + (V>,0>) ^ 



' l Je (s)0(x)X £ (.s,a;) J 'dW A (s,a;), 



i) 



X £ G C(R+,C^ 
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Similar reasoning gives 



(3.12) 



{Y t s A) = {A%{t), 



a; 



A 
~2 



n (f>) + (i/>, 4>) ds 



l Je (s)<l>(x)Y £ (s,xrdW( S ,x), 



Since J7J 



that 
(3.13) 



r eC(R + ,C+ p ), 

- T e +e, the alternating summation in the definition of implies 



sup \A E (t,x)\ < etp(x) 



It follows from Q3.2p and (J376J) (recall that 6 = f ip(x)dx) that 

X e (t,a;)ljc(t) < 2 / P s 4>(x)ds < Ae 1/2 b. 
Jo 

Therefore for any T > and <f> as above 



(3.14) 



\ja{s){X e {s,x)Ycj){x)dW{s,x) ) < {4e 1/2 b) 2p T\ 



By identifying the white noise W with associated Brownian sheet, we may 
view W as a stochastic process with sample paths in C(R+, Ct em (]R)). Using 
bounds in Section 6 of |Shi94j (see especially the pth moment bounds in the 
proofs of Theorems 2.2 and 2.5 there) it is straightforward to verify that for 
e n i 0, {{X^,Y^,W) : n e N} is tight in C(R+, (C+ p ) 2 x C tem ). Some of 
the required bounds are in fact derived in the proof of Lemma [6] below. By 
(I3.13|) , (|3.14[) and their analogues for Y £ , one sees from (|3.1ip and (|3.12[) that 
for any limit point (X, Y, W) , X and Y are C+j p -valued solutions of (|3.1j) with 
respect to the common W. It remains to show that X and Y are distinct. 

We know X e (t, •) and Y e (t, ■) will be locally Holder continuous of index 
1/4 but it will be convenient to have a slightly stronger statement. We note 
parenthetically that any other index of Holder continuity for X s (t, •) and 
Y e (t, ■) would yield the same range for p in Theorem provided that the 
index were less than 1/2. Let 

= inf{s > : 3x, x G R such that 

\X £ (s,x) -X e (s,x')\ + \Y e (s,x) -Y e (s,x')\ > k\x - x'| 1/4 }. 



We will show in Lemma [6] that lim; c _j. 00 
M 6 N. 

Note that on [T£j, U^), Y e = and 
\X £ (t,x)-X £ (t,x')\ = 2 



sup 0<£<1 P(V£ < M) = for any 



Pt - T e.(z){ip(z+x)-^(z+x'))dz < 2\\i>'\\ 00 \x-x'\. 
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This implies that on the above interval for all real x,x' , 



\X s (t, x) ~ X%t, x')\ + \Y e (t, x) - Y s (t, x')\ < 4(|| V'Hoo V U\\oo)\x - xf/ 4 , 



where the inequality holds trivially for \x — x'\ > 1 since the left side is at 
most 4IIVH oo . By symmetry it also holds on Plj+n We may assume 

k > 4(||-0'|| oo V 1 1 V 1 1 oo ) and so the above implies 



(3.i5) v k E € U(c/;,t; +1 ]u{oo}. 

3=0 



We fix a value of k which will be chosen sufficiently large below. We will now 
enlarge our probability space to include a pair of processes (Xf, Y t 6 ) which will 
equal {{X e , 1), (Y t £ , 1)) up to time V£ and then switch to a pair of approximate 
solutions to a convenient SDE. Set p' = ^ € (0, \ ) and K(k) = E^l/A, k) 
where -ftjg] is as in Theorem [5J We may assume our (D,,J-,J- t ,P) carries a 
standard ^-Brownian motion (B s : s < V£), independent of (X s , Y £ , W). 

To define a law Q a on C(R, M+)x [0, oo], first construct a solution (Y e , D e ,B) 

of 



(3.16) 

Y t s = b [ l(e <s)ds+ [ l(e < s)(Y s e f V ' K{k) dB s , Y e > 0, 
Jo Jo 

Dl=2b(tAe)+ f\(e < s)[(Y: + Dlf' - (Y^r']^K{k)dB s , D s > 0. 
Jo 



Such a weak solution may again be found by approximation by solutions of 
Lipschitz SDE's as in Theorems 2.5 and 2.6 of |Shi94j for the more complicated 
stochastic pde setting. Set X e = Y £ + D £ and ff = inf{i : X e t = 0} > e, and 
define 



(3.17) 



Qo{A) = P{{X^ fr Y^ fr Tt)EA). 
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Next we enlarge our space to include (X £ ,Y £ ) so that for finite t < Tf 
(this time is defined below), 

(3.18) Yf = {Y& Vj .,l) + b(t-V£) + 

+ f l(s>Vi)(Y s e ) p '^W)dB s , Y £ >0, 
Jo 

(3.19) D$ = (DI AV s,l) 

+ [ l(s > V k £ )[(D e s + W - (Y s e ) P ']VW)dB s , D £ > 0, 
Jo 

(3.20) XI = Y t * +Df = (X? AV * , 1) + b(t - V£) + 

+ [ l(s>V£)(Xiy'^K(k)dB s , 
Jo 

(3.21) Tf = inf{i > : X? = 0} < oo. 

Note that if V £ > Tf, then X e = (X £ ,l) for_t < Tf and so Tf = Tf. 
Therefore, Tf A V£ < Tf . We conclude that Xf AV e = Y £ AV , + D £ AV e for 
t < Tf , thus proving (f3~20]) . 

To carry out the above construction first build (Y V e +t , Dy^ +t , By^+t — 
Bv e ) by approximation by solutions to SDE's with Lipschitz coefficients as 
in Theorem 2.5 of |Shi94] . This and a measurable selection argument (see 
Section 12.2 of jSVj ) allows us to build the appropriate regular conditional 
probability 

Q°(Y^,l),(D' v ,,l){-) 

ee P(((Y v , +t: DU +t ,B(V£ + i) - B(VZ)),t > 0) G -\x £ , y £ , w), 



x 



. V B +t ,U V B +t , 

where {Q® d ■ y, d > 0} is a measurable family of laws on C(R+,R+ 
This then allows us to construct (X s ,Y £ , D £ ) as above on an enlargement 
of our original space which we still denote (O, J 7 , J-"t,P)- We also may now 
prescribe another measurable family of laws {Q y . x ■ (y,x) € C(R+,1R+)} on 
C(R+,K+) x [0, oo] such that for each Borel A, w.p. 1, 

(3-22) Q {x . ■ 1}i{Y . 1) {A)=P{{X* At Y Af T[)eA\X*,Yt,W). 

k k 

Define 

(3.23) Qi(A) = P((X': aT Y Af Tf) e A) = e(q {x . 1)t{Y% 1} (A) 

1 1 \ ,xv k k 

Next, inductively define (Xf,Y t e ), t g [Tf,Tf +1 ], and {Tf} in a manner 
reminiscent of that for (X £ ,Y £ ), and consistent with the above construction 
for j = (set Tg = 0). Assume the construction up to Ty is such that 

(3.24) X £ Tl = Yf, = on {Tf,. < oo}, 
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and 

(3.25) % = Tl } on {V§ > %}. 
Define 

(3.26) Xj (t) = ^(y;+t)AT« +1 , 
and similarly define Yj. On {TJ^ < V£} set 

(3.27) P((X 2j ,Y 2j ,fl j+1 - flj) G \T n . V a(X°, Y s , W)) 

On {oo > f| 3 - > V£} set 

(3.28) P((X 2j) y 2j -,f| i+1 - %) G V cr(X e , Y^ W)) = Q o (0- 

These definitions imply TJ^x = inf{t > 2% : X t e = 0} and that on our 
enlarged probability space, conditional on Jy e and on {V£ > 2^}, (|3.18l) - 

®M hold for * € while on {oo >% > V£}, for t G [% T| j+1 ], 

(I3A61) . (I3~T71) and (13^251) give 

Y t E = b [ 1(2% +e< s)ds 
Jo 

+ / l(2f . + e < S )(F// ^K{kjdB( S ), Y E > 0, 
Jo 

(3.29) 

D° t =2b{(t-T$ j )Ae) 

+ f 1(2% + e < S )((F/ + 5^' - (W) V 7 ^) dB., & > 0, 
Jo 

X t £ = 2b((t - 2%) A e) + b [ 1(2% +e<s)ds 

Jo 

1(2% + e < s){Xtf' ^/K{k)dB{ s ). 



Now assume 2% +1 < oo and construct (X^,Y t £ ) and D\ = \X^ — Y t £ \ for 
^ e Psy+i> -^2.7+2] as a bovc but with the roles of X and Y reversed. This 
means that on {2% +1 < V£}, 

(3.30) P((Y 2j+1 ,X 2j+1 ,fi j+2 - 2% +1 ) e -| V ff (I £ ,r, W)) 

= g ^i J+1+ -)-r 1> -^H +1 +-)-r 1)(,) ' 

and on {00 > 2% +1 > V fc e }, the above conditional probability is again Qq. 
The apparent lack of symmetry in the definitions arises because we have also 
reversed the roles of X e and Y e on [2% +1 , 2% +2 ] . The above definition implies 
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that Df = Yf -Xf>0 on [% +1 ,% +2 ], T| J+2 = inf{t > Tf j+1 : F t £ = 0}, 
and X E (T|. +2 ) = F £ (7% +2 ) = 011 {^+2 < °°}- 

It follows from (pTT8]) . (pT20)) (now with i G [7%, T| j+1 ]) and (|3~25j) that on 
{Vu > we have f. 



T, 



2j + l- 



2j+1 — - L 2j+i- Symmetric reasoning shows that on 
i v k > Ti j+2 }, Ti j+2 = Ti j+2 . We have verified (j3"^4]) and for j + 1. 

Since f/ +1 - f/ > e (by . (f3T20|) and (|3^9)» . f| t °o and our inductive 

definition is complete. 

The reasoning above to show T^j^i 
obvious induction also shows that 



Tfi+i on {V£ > T|, +1 } and the 



(3.31) 



AW = I (-^tAVi . !> - O'tAV? ,1)1 Vt > a.s. 



The following consequence of the above construction will be important for 



P((X 2 j,Y 2 j,T2 j+1 - T 2j ) G - 



(3.32) 



and 



(3.33) 



P((Y 2j+ i,X 2jA 



h(-) a.s. on {5^<^f} 
T 2 j+i ) G ■|j 7 T't 



2j+2 



h(-) a.s. on {T|. +1 <F fe £ } 



= Qo(-) a.s. on {F fe e < 2%}, 
P((?2j+l)^2j+i,r| j+2 - T| j+1 ) G •|J 7 f« f+1 ) 
= Q (.) a.s. on {V£ < T S J+1 }. 

Consider, for example, the first equality in (|3.32[) . By p.27|) we have for a 
Borel set B and A G T^*. , A c {Vjf > 2%}, 



P({(X 2j , Y 2j ,Tl j+1 ~ %) G B} n A) 



B<3( 



\ 



„i>( B )L 



(3.34) 



(B) 



Ft- 



AV' 



In the last line we used the fact that VJ: > T 2 j 



T. 



on A to see that 



A G Ft*.kV£- Formula ([3T3T|) shows that our construction of (X s , Y £ ) has 



not increased the information in J^t^.av so we may use (|3 . 9[) . Applying Q3.9p 
and the fact that = T^- + V£ o T e on {V" fc e > T^}, where (0 t ) arc the shift 
operators for (X s , Y s ), we conclude from (|3.34p that the far left-hand side of 
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(|3~34| equals 



E 



(Q(x° Ave ,i>,(i? AV . ,i> ( B )) U) = B (Qi , 



by (|3.23p . This gives the first equality in (|3.32p and the second inequality 
holds by a symmetric argument. The proof of (|3.33[> is easier. 

Our next goal is to show there is positive probability, independent of e, of 
D £ hitting some appropriately chosen x G (0, 1) before X s or Y £ hits 1. By 
(|3.32p and (|3.33p the excursions of X s V Y £ away from are governed by Qo 
or Q±, depending on whether or not V£ has occurred. Therefore we need to 
analyze these two laws. 

Consider the more complex Q\ first. Use (|3.1ip . with tj> = 1, in (|3.20p and 
the fact that V% > e (by jf3~T5]) ) to conclude that under Q 1: X £ = 2b(t A e) 
for t < e and for < t < ff - e, 

X; +e =eb+ ((t + e) A V§)b + b f + * l(s > V^) ds 

t+e 



I) 



t+e 



(3.35) 



l(s < Vg)X e (s,x) p dW{s,x) 

+ y ' l( S >^ £ )(X s T'v^(fc)^ s 
2eb + tb + N t , 



where N is a continuous {Tt+e)— local martingale such that 



1(* < ^ £ ) / ^(a.a)^ + l(s > Vg)(Xtf>'K{k) 



ds 



t+e 



{N)'(s)ds. 



By the definition of V fc £ we may apply Theorem [5] with (a, j3, C) = (2p, 1/4, k) 
and conclude that 



(3.36) (N)'(s) > l(s<V£)K{k)U X £ (s,x)dx 

+ l(s>V k s )(X £ s fP'K(k) 
= K{k){Xtfv', 



((p/2) + l)/(5/4) 



where (|3.3ip is used in the last line. 
Define a random time change r< by 



(3.37) 



t = 



(NY(s + s) 

o K{k){x £ s+e yp 



7 ds = A(T t ), t<A(T£-e). 
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The restriction on t ensures we are not dividing by zero in the above integrand 
because T-f is the hitting time of by X s . Clearly (|3.36[) implies 

(3.38) r'(t) < 1 for t < A(T[ - e). 
For t < A(Tf -e), let 

(3.39) X(t) = X%T t +e) = 2be + br(t) + N(t), 
where N t = N(r t ) is continuous (J r n+e )-local martingale such that 

(N)t= f * + (N)'{s)ds= I 1 (NY s+£ ds = K(k) f X(r) 2p ' dr. 

Js JO JO 

This follows by using the substitution s — r r and calculating the differential 
dr(r) from (pOT)) . Note also that if f x = ini{t > : X(t) = x}, then 
A(T[ — e) = Tq. Therefore by (|3.39[) we may assume there is a Brownian 
motion B so that 

(3.40) X(t A To) = 2fe + br(t A T ) + / yj K{k)X(s) p ' dB(s). 



The scale function for a diffusion defined by a similar formula, but with t A To 
in place of r(t A To), is 



/" 2by 1 - 2 P' I , 

Sk{x)= J CXP \-K(k)(l-2pi)\ dy - 

That is, Sfe satisfies 

K(k)x 2p ' 



(3.41) 4(x)+bs' k (x) = on [0,oo), s k (0) = 0. 

By Ito's Lemma 

f tAfo - K(k)X(u) 2p ' 

s k (X(t A T )) = s k (2be) + J br> ' (u)s' k (X(u)) + [ > ^ 1 4{X{u)) du 

tf\f a 

s' k (X(u))dN(u). 

p.38p and (|3.41[) show that the integrand in the drift term above is non- 
positive, and so Sk(X(tAToATi)) is a supermartingale which therefore satisfies 

E(s k (X(tAf AT 1 )) < s k {2be). 

This implies that 

Qi{X e t = 1 for some t < Tf) = Qi(s k (X(- A f x A T )) hits s k (l) before 0) 

= lim QJskiXitAfoAT^/s^l)) 

(3.42) < Sk (2be)/s k (l). 
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Under Qq add the equations in p,16|) to see that (we write X s for X e ), 

X E t+£ = 2be + bt + J {X E S Y yjK{k) dB s , t + s<T[ = inf{f : Xf = 0}. 

This is equation (|3.40j) with t in place of r t and so the previous calculation 
applies to again give us (|3.42[) with Qq in place of Q\. 
Under either Qi, Xf = XI V Y t e and so we conclude 

(3.43) Q % {X e t V Y t £ hits 1 for t < Tf ) < s k (2be)/ s k (l), i = 1,2. 

We next consider the escape probability for D e under Q\, Let xq G (2be, 1) 
and 

T D (Q,x Q ) = inf{< : Df = or x } < Tf Qi - a.s., 
the last since D\ = X e t - Y t e < X e t for t < Tf under Q x . It follows from 
(gH) and (j3~T§)) that £)f = 2b(t A e) for t < e, and for £ + e < Tf we have, 

(3.44) r^+e = 26e + / f (X%s, xf - Y s {s, < V£) x) 



+ ^ t+£ i( s > ^)((X 8 T' - {Y;Y')^K{k)dB s , 

which is a non-negative local martingale in t. We have 
Qi(3i<Tf: A e >^o) 

> T^(T 5 (0, 2 ; o )) a ;- 1 l(T B (0,x o ) < oo.Tf < oo)) 

(3.45) = e(D s (T d (0,x ) A 7\ E )io ') - s(5 £ (T 5 (0,x ))xo ^(Tf = oo)). 

The first term on the right-hand side is the terminal element of a bounded 
martingale and so 

(3.46) E(D s (T D (0,x Q )Af^x^ = 2bs/x . 

It follows from (|3~3"5j) that on {Tf = oo}, 

(3.47) {limsupXf < oo} C { lim N t = -oo}, 

t-yoo t^oc 

which is a Qi-null set by the Dubins-Schwarz theorem which asserts that a 
continuous martingale is a time-changed Brownian motion. Therefore 

(3.48) X s t+e hits or 1 for some t < Tf - e, t < oo, Qi - a.s., 
and therefore 

(3.49) e(D s (T b {0,x ))xo H(Tf = oo)) < Qi(Tf = oo) 

< Qi(X% hits 1 for t < Tf) 

(3.50) < s fc (2fe)/s fe (l), 
the last by ([g^]) . 
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Since lim x _>.o+ Sk{x)/x = 1, there is an eo(fc) > and xq = xo(k) <G (0, 1), 
such that 

(3.51) e < £o implies Sk(2be) < 3be and 26e < xo < Sfc(l)/6. 

So for £ < £o an d xo as above we may use (|3.46j) and (|3.49[) in (|3.45p and 
conclude 

Qi(3t G [£,ff) : S? > x ) > (2bs/x ) ~ (a fc (2te)/« fc (l)) > (&e)/x . 

Virtually the same proof (it is actually simpler) works for Q$. Under Qi, 
D e t = \Xf - Y t e \ for t<ff = inf{t : X e t V F t e = 1} and so we have proved for 
xq as above, 

(3.52) Q l (3t G [e, ff] : \XI - F/| > zo) > — for i = 1, 2 and < £ < £ , 

Xo 

and (see (pT48| for i = l) 

(3.53) X t £ V hits or 1 for t < Tf, t finite Qi - a.s., i = 1, 2. 
Let 

M = min{j : (X e V F £ )(t + T/) hits 1 for f < f/ +1 - T/}, 

and 

Mi = min{j : \X e - Y E \(t + ff) hits x for t < T/ +1 - Tf}. 
Use (f3T32|) . (pT33l) and (|3~43f to sec that 
P(M > n) 

= > n - 1)P(X £ V y e ((^_i + •) A T £ ) doesn't hit 1^?*^ 

s k (2b, 



> P(M >„-i)(i-fig2) 



Therefore, if p x = £ ff|fj i , then 

(3.54) P(A/i > n) > (1 -pi) n+1 . 

Similar reasoning using (|3.52[) in place of (|3.43|) shows that if p<i = then 

for £ < £o, 

(3.55) P{M 2 >n) < (1- P2 ) n+1 . 
Note that (|3.51[) shows that 

(3.56) E1 = ^^ S J^1> 1 ^1> 2 . 
pi s k (2be) x 3 x 



18 K. BURDZY, C. MUELLER, AND E.A. PERKINS 

If n = [p^ 1 ] we get for e < eq 

P(JV 2 < M) > P(M > n) - P{M 2 > n) > (1 - Pl ) n+1 - (1 - p 2 )" +1 

>(l-p 1 )"+ 1 -(l-2p 1 )"+ 1 

^(e^-e" 2 ), 

where the last inequality holds by decreasing £o(k), if necessary. If 

t £ = inf{i : X e t V Y t E > 1}, 
then the above bound implies that for £ < £o, 

(3.57) P(sup \XI - Y t °\ > so) > ^(e- 1 - e" 2 ). 

t<t e 1 

Now let 

t £ =inf{i:(X|,l}v(y/,l}>l}. 
Then (f33Tj) shows that 

if t 6 < V£, then t 6 = t e and (X { £ , F/) = {{XI, 1), (F/, 1)) for all t < t e , 
and so by (|3.57j) for e < £ , 

(3.58) P(sup |(X t £ , 1) - (Y/, 1)| > x ) > ^{e- 1 - e" 2 ) - P(V k E < t e ). 

Now recall we have e„ | so that (X £ ™ , F e " , W) -> (X, F, VP") weakly on 
C(R+, (C+j p ) 2 x Ctem)) where X and F are C^ p -valued solutions of (|3.1|) . 
Arguing as in (|3.47|) and using Dubins-Schwarz, we see that 

(3.59) limsup(X t , 1) = limsup(F t , 1) = oo a.s. 

t— >oc t— foo 

Standard weak convergence arguments now show that {t En } are stochastically 
bounded. Lemma |6] therefore shows that we may choose a fixed k sufficiently 
large so that 

P(V k En < t e J < -^(e" 1 - e" 2 ) for all n. 
Using this fixed k throughout we see from (|3.58[) that for large enough n 

P^sup \{X £ ",1) - (F/M>| > x ) > \{e-' - e~ 2 ). 

If t' = in£{t : {X t , 1) V (F t , 1) > 2} < oo a.s., by ([339]) . then the above implies 

p(sup|(X t) l)-(F t ,l)| >x /2) > -(e^-e- 2 ), 
H<t' ' 4 

and so P{X ^ F) > ±(e _1 - e~ 2 ). □ 
Lemma 6. For any M e N. limfc_ ) . 00 sup 0<e<1 P{V k e < M) = 0. 
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Proof. The proof depends on a standard argument in the spirit of Kolmogorov's 
continuity lemma, so we will omit some details. 
Fix the time interval [0,M]. Define 

Vg(X) = inf{s > : 3x, x' G K such that 

\X E {s,x)-X E {s,x')\ > k\x-x'\ 1/4 }. 

and V£(Y) likewise. It suffices to prove Lemma [6] for V fc e replaced by V£(X) 
and V£(Y) and so clearly we only need consider V£(X). Recall that {X e ™} 
is tight in C(R+, C+ ap ). So it suffices to choose a constant K > and prove 
the lemma for X 6 (t, x) A (Ke~\ x \) 1 / p in place of X £ . Considering the integral 
equation for X s , and using the fact that ip G C^(R) we see that it is enough 
to prove Lemma [5] with X 6 replaced by the stochastic convolution 

N e {t,x) = [ p t -s(x - y)ip £ (s,y)dW(s,y). 



Here one can use Lemma 6.2 of |Shi94] to handle the drift terms. The term 
f e (s, y) is a predictable random field satisfying 

\<p e (t,x)\ < Ke~W 

for all t € [0, M], x G R almost surely. Since our estimates are uniform in e, 
we will omit the superscript on tp and N from now on. The constants below 
may depend on M and K. 

Now we rely on some standard estimates which are easy to verify. We claim 
that there exist constants qo, Kq such that for0<t<£ + (5<-^ and ieR. 
and for S < 1, 

(3.60) / / p 2 s (x - y)e- 2pM dyds < RaSh'^, 



JO JR 

{pt-s+s(x-y)-Pt-s(x-y)} 2 e- 2pM dyd S < K 8ie- q °W, 

JR 

\pt- s (x-y + S)-p t . s (x-y)] 2 e- 2 ^dyds < K 6e~^. 

10 JR 

From these inequalities, it follows in a standard way that for some positive 
constants q\ , Co , C\ , we have 

(3.61) P(\N(t + 6,x) — N(t,x)\> X) < C exp(-^-\ e -« l|x| , 

P{\N(t,x)-N(t,x + 5)\>\) < C exp(-^^\ e ~^. 

For example, if we write 

M r = / / [pt- s (x~y + 5) - p t - s (x - y)]ip(s, y)W(dy, ds) 
Jo JR 
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then M r is a continuous martingale and hence a time changed Brownian mo- 
tion, with time scale 

E(r) = [p t ~s(x - y + 5) - p t -s(x - y)] 2 tp 2 (s,y)dyds 

JO JR 



< 



f / \p t _ s (x-y + 6)- Pt _ s (x-y)} 2 K 2 e- 2 Mdyds. 

Jo JR 



Thus, 

P(\N(t,x) — N(t,x + 5)\ > A) = P(|M t |>A) 

< P\ Sup |fl a |>Aj 

\0<s<E(t) J 

and then the reflection principle for Brownian motion and the third inequality 
in (|3.60p (to bound E(t) for t < M) gives the second inequality in (|3.61|) . 

Now we outline a standard chaining argument, and for simplicity assume 
that M = 1. Let Q n be the grid of points 



k t 

22n ' 2" 



< k < 2 2n ,£e Z 



The Borel-Cantelli lemma along with (|3.61j) now implies that for large enough 
(random) K\, if n > K\ and p\,p2 are neighboring grid points in Q n , then 

(3.62) |JV(pi)-JV(pa)|<2-*. 

Now suppose that = (i,-, Xj) with (xj — x 2 | < 1, and that each point % lies in 
some grid C?„ . From the above, there is a path from q\ to (72 utilizing edges in 
grids Q n i, with n' < n, each edge in the path being a nearest neighbor edge in 
Q n ', and with at most 8 edges from a given grid index n' . Let no be the least 
grid index used in this path. We claim that for some constants C > c > 0, 
such a path exists with no satisfying 

c2 -2n < | tl _ t2 | < C2" 2 ™ , 

c2- no < \ Xl -x 2 \ < C2~ no . 

Using the triangle inequality to sum differences of N(t, x) over edges of the 
path, we arrive at a geometric series, and conclude that 

(3.63) \N( qi ) - N(q 2 )\ < Ci2~^ if n Q > K x . 

Although we have only proved the above for grid points, such points are 
dense in [0,T] x R, and N(t,x) has a continuous version because X(t,x) is 
continuous, and the drift contribution is smooth. Therefore it follows for all 
points in [0, 1] x K. We have proved ([5jB5} for ||g x - q 2 \\ < C2" A ' 1 where K x 
is stochastically bounded uniformly in e. The required result follows. □ 



NON-UNIQUENESS FOR SPDE 



21 



Sketch of Proof of Theorem [4l We carry out an excursion construction of 
an approximate solution X s to (|1.2j) by starting the ith excursion at (— lYetp, 
and then run each independent excursion according to a fixed law of a C^ ap - 
valued solution to (jl.2p with X = etp, if i is even, and its negative if i is odd, 
until the total mass hits 0. At this point a new excursion is started in the 
same manner. Theorem [5] is used to time change X^{1) into an approximate 
solution Y £ (t) = X^e(l) of Girsanov's equation 

(3.64) dY t = \Y t \P'dB t , 

with p < p' < 1/2 and dT J < 1- There will be an additional term A E (t) 
arising from all the excursion signed initial values up to time t but it will 
converge to uniformly in t due to the alternating nature of the sum. We 
now proceed as in the excursion-based construction of non-zero solutions to 
Girsanov's sde (|3.64[) to show that one of the excursions of the approximate 
solutions will hit ±1 before time T with probability close to 1 as T gets large, 
uniformly in e. Let N e be the number of excursions of Y e until one hits ±1 
and let N e (T) be the number of excursions of Y s completed by time T. N e 
is geometric with mean e _1 by optional stopping. Let Ui(e) be the time to 
completion of the ith excursion of Y £ . Assuming \fTe~ 1 € N, wc have 

P(sup|F/|>l)>P(iV e (T)>iV £ ) 

s<T 

> P{N e {T) > VTe- 1 ) - P(N £ > VTe- 1 ) 

(3.65) > P(C/^ e - 1 (e) <T)- {l- e )^\ 

A key step now is to use diffusion theory to show that if Y satisfies (|3.64p 
(pathwise unique until it hits zero) then 

(3.66) P(Y t > for all t < T\Y = 1)) - c T~ 1/(2(1 - p)) as T -> oo. 

If Ui (1) is the time of completion of the ith excursion of Y where the excursions 
now start at ±1, then scaling shows that 

P{U^ 00 < T) = PiU^ (1) < e-^T). 

(|3~66| shows that f/^-i (l)/(\Zre" 1 ) 2 ( 1 - p ) converges weakly as e | to a 
stable subordinator of index a = (2(1 — p))^ 1 and so for any rj > we may 
choose T large enough so that for small enough e (by (|3.65p ) we have 

P(sup |F/| > 1) > P(U^ (e) < T) - (1 - e)^ 1 

s<T 

>p ( Uyre-^) <TP ]_ e -Vr >1 _ 

The fact that {T £ )'(t) < 1 allows us to conclude that with probability at least 
l—i], uniformly in s, the total mass of our approximate solution X s will hit 
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±1 for some t < T. By taking a weak limit point of the X 6 we obtain the 
required non-zero solution to (|1.2j) . □ 

4. Pathwise Non-uniqueness and Uniqueness in Law for an SDE. 
The stochastic differential equation corresponding to (|3.1[) would be 

(4.1) X t = X Q + bt + f {X s ) p dB s , X t >0Vt>0 a.s. 

Here 6 > 0, < p < 1/2, B is a standard (J^-Brownian motion on (O, T , T t ,P) 
and X is J-" -measurable. A much simpler argument than that used to 
prove pathwise non-uniqueness in (|3.1[) allows one to establish pathwise non- 
uniqueness in (|4.1[) . One only needs to apply the idea behind construction 
of (Xt,Y t e ) for t > V£. In any case the result is undoubtedly known, given 
the well-known Girsanov examples (see, e.g., Section V.26 in |RWj ) and so we 
omit the proof. 

Theorem 7. There is a filtered probability space (£l,J r ,J 7 tl P) carrying a 
standard T t -Brownian motion and two solutions, X 1 and X 2 , to (|4.1[) with 
Xl = Xl = X = such that PiX 1 ^ X 2 ) > 0. 

Weak existence of solutions to (|4.1[) for a given initial law may be con- 
structed through approximation by Lipschitz coefficients. This is in fact how 
Shiga |Shi94j constructed solutions to (1.1) and hence is the method used in 
Theorem 1. As we were not able to verify whether or not uniqueness in law 
holds in (|3.1j) it is perhaps interesting that it does hold in (|4.1j) . That is, the 
law of X is uniquely determined by the law of Xo . We have not been able to 
find this result in the literature and since the solutions to (|4.ip turn out to 
have a particular sticky boundary condition at which was not immediately 
obvious to us, we include the elementary proof here. 

Theorem 8. Any solution to (14. ip is the diffusion on [0,oo) with scale func- 
tion 

(4.2) , W = £e*p{;S^> 
(with inverse function s _1 on [0,s(oo)J, speed measure 

(4.3) m(dx) = ^-ifc^-i^p + b^Soidx) on [0, s(oo)), 

and starting with the law of Xq. In particular if Tq = ini{t : X t = 0}, then 

/•To+e 

(4.4) T < oo implies / 1(X S = 0) ds > Ve > a.s., 

J T 

and solutions to (|4.1[) are unique in law. 
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Proof. The last statement is immediate from the first assertion. 

To prove X is the diffusion described above, by conditioning on Xq we may 
assume Xq = xq is constant. We will show directly that X is the appropriate 
scale and time change of a reflecting Brownian motion. Note that s is strictly 
increasing on [0, 00) (in fact on the entire real line) so that s _1 is well-defined. 
Note also that 

r -2b\x\ 1 ~ 2p ] 

s'(x) = exp| — J is of bounded variation and continuous, 

and 



(4.5) s"(x) 



-s'(x)2bx- 2 P ifx>0, 
a' (x)2b\x\~ 2p ifx<0. 



If Lf(x) is the semimartingale local time of X, Meyer's generalized ltd formula 

(see Section IV.45 of |RWj ) shows that 

(4.6) 



Y t = s(X t ) =Y + / s'(X u )X p dB u + b s\X u ) du+- Lf {x)ds'{x). 



1 



Since s' is continuous at 0, 

(4.7) \ J Lf{x)ds'{x) = \Jl{x> 0)L?(x)ds'(x) 

= \\ 1( - X> °) L ?( x ) s "( x ) dx 

= -b [ s'(X u )X- 2p Xl P l(X u > 0)du (by 63D) 
Jo 

= -b f s'(X u )l(X u >0)du. 
Jo 

So (H) and s'(0) = 1 imply 

(4.8) Y t = Y + f s'(X u )X p dB u + b f l(X u - 0) du. 

Jo Jo 

Define U = J °° s' '{X u ) 2 X 2p du and a random time change a : [0, U) — > [0, 00) 

by 

ra(t) 

(4.9) / s'(X u ) 2 X 2p du = t. 

Jo 

Clearly a is strictly increasing and is also continuous since X cannot be on 
any interval. If R(t) = Y(a(t)) for t < U we now show that R is a reflecting 
Brownian motion on [0, 00), starting at Yq, where we extend the definition for 
t > U by appending a conditionally independent reflecting Brownian motion 
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starting at the appropriate point. In what follows we may assume t < U as 
the values of R(t) for t > U will not be relevant. We have from (|4.8j) 

,a(i) 

R{t) = p t + b l(X u = 0)du = f3 t + A u 



where (/3) t = t and so /3 is a Brownian motion starting at Yq. A is contin- 
uous non-decreasing and supported by {t : X(a(t)) = 0} = {t : R(t) = 0}. 
By uniqueness of the Skorokhod problem (see Section V.6 in [RWj ) R is a 
reflecting Brownian motion and A t = (0), that is, 

/•«(*) 

(4.10) b / 1{X U = 0)du = Lf (0). 



Let a 1 : [0, oo) — > [0, U) denote the inverse function to a. Now differentiate 
P~!?)) to see that 

(4.11) if X u > 0, then (aT 1 )'^) = s'(X„) 2 X^ . 

We may use (|4.10|) and (|4.11|) to conclude that 

t= f l(X u >0)du+ [ l{X u = 0)du 
Jo Jo 

= fo d{a ' 1{u)) + b ~ lL *-Ht)(°) 

r-- X W l(R(v)>0) R 

s>(s-i(R(vWs-i(R(v)yp d +b L "-H*)^' 
where we have set u = a(v) in the last. Therefore if m is as in (|4.3[) . then 

t = / L R ( a - 1 (t),x)dm(x), 

J[0,cx>) 

and 

X(t) = s- 1 (R(a-\t)). 
This identifies X as the diffusion on [0, oo) with the given scale function and 
speed measure. □ 

Remarks. (1) One can also argue in the opposite direction. That is, given 
a diffusion X with speed measure and scale function as above and a given 
initial law on [0,oo), one can build a Brownian motion B, perhaps on an 
enlarged probability space, so that X satisfies (|4.1[) . giving us an alternative 
weak existence proof. 

(2) One can construct solutions as weak limits of difference equations or cquiv- 
alently as standard parts of an infinitesimal difference equation. Here one cuts 
off the martingale part when the solution overshoots into the negative half- 
line and lets the positive drift with slope b bring it back to R+. The smaller 
the b the longer it takes to become positive, the more time the solution will 
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spend at zero and so the larger the atom of the speed measure at 0. A short 
calculation shows that at p = 1/2 the overshoot reduces to At (the time step 
in the difference equation) and so there is no time spent at in the limit. (See 
Section V.48 of |RWj for the standard analysis.) 

(3) It would appear that (|4.1|) is not a particular effective tool to study diffu- 
sions with drift b on the positive half-line. By just extending the equation to 
[0, oo ) we inadvertently pick out a particular case of Feller's possible boundary 
behaviors at among all diffusions satisfying (|4.1[) on (0, oo). (This is cer- 
tainly not a novel observation — see the comments in Section V.48 in }RWj .) 
Presumably things can only get worse for the stochastic pdc (|3.1[) . In the next 
section we scratch the surface of this issue and show that all solutions to this 
stochastic pde spend positive time in the (infinite-dimensional) zero state. 

5. Proof of Theorem [3] 
Let X be a solution of (|3.1|) and define 

V k = inf{s : 3a/, x such that \X(s,x) - X(s,x')\ > k\x - x'\ 1/4 }. 

As in Lemma [6] (but as there is no e it is a bit easier), lim^ Vk = oo a.s. As 
in Section El we set p 1 = ^ and K(k) = i^l/4, k). If we define 

R[u) = J Amxt)\%' ^(X(u),l)>0,ndu<V k , 
[ 1 otherwise, 

then by Theorem [5l 

R(u) > 1 for all u > 0. 
Introduce a random time change r given by 

r(t) 



(5.1) / R(u)l({X u ,l) >0) + l((X u ,l) 



0) du = t. 



Clearly r is strictly increasing, continuous and well-defined for all t > 0. 
Differentiate (|5.ip to see that 

T'(t)R T(t) l((X T(th 1) > 0) + r'(i)l((A T(t) , 1) = 0) = 1 for a.a. t > 

(a. a. is with respect to Lebesgue measure), and therefore 

(5.2) r'(t) = R;f t) l((X Tit) , 1) > 0) + l((A T(t) , 1) = 0) < 1 for a.a. t > 0. 

Now let 

Y(t) = (X(r(t)),l)=br(t) + M(t), 
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where M is a continuous local martingale satisfying 
r-r(t) 



(M) t = J X(u,x) 2p dxdu 

X(T{r),x) 2p dx\R{T{r))- 1 l{Y(r) > 0) + l(Y(r) = 0) 



dr 

K(k)Y 2p ' dr for r(t) < VJfc. 

We have used (|5.2| in the second line. If Tfc = r _1 (Vfc) (a stopping time w.r.t 
the time-changed filtration), we may therefore assume there is a Brownian 
motion B so that 

y(i A T fc ) = fcr(t A T fc ) + / ^K(k)Y p dB r . 



If b' = Kikf/W-^b, then F(t) = if ( fc ) 1/(2(p '~^Y{t) satisfies 



F(f A T fc ) = &V(t A T fc ) + / 1? dB 



tAT fc 



If s(x) = cxp| 2b i^ 2p , — | dy, then an application of Meyer's generalized 
Ito's formula shows that if Z(t) = s(Y(t)), then for t < Tj., 

Z(t) = J s'(Y r )Yf dB r + V j s'{Y r )r'{r) dr + ^J Lf(x)ds'(x). 

Here, as before, L Y is the semimartingale local time of Y. Now argue as in 
dHZH to see that for t < T k , 



i J L?(x)ds'(x) = -b' J s'(Y r )l(Y r > 0)dr. 

Therefore if N(t) = J* s' (Y r )Y p ' dB r and 

A{t) = b' [ s'(Y r )(l - r'(r))l(5> > 0) dr, 
Jo 

then for t < T k , 

Z{t) = N{t) - A{t) + b' [ s'(0)r'(r)l(r(r) = 0) dr 
Jo 

(5.3) =N{t)-A(t) + b' j l(Y(r)=0)dr, 



where we used s'(0) = 1 and (|5.2p in the last line. A is a non-decreasing 
continuous process by (|5.2[) . N is a continuous local martingale, and N(0) = 
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A(0) = 0. Fix k and assume V% > 0, and so T k > because T k > Vk- If 



then by (|5.3j) . Z(t A T+) is a continuous non-negative local supermartingale 
starting at and so is identically zero. This means Z(r) = for r < T + and 
so the same holds for Y(r), which by the definition of 7+ and assumption that 
Tj. > implies that T + = a.s. Since t °o a.s. we have shown that w.p. 1 



Setting r(r) = u and using (|5.2p again (to show r'(r) = 1 on {(X(r(r)), 1) = 
0} for a. a. r), we see that the above implies 
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